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ABSTRACT 



Unravelling the morphologies of Luminous Compact Galaxies 

S ' using the HST/ACS GOODS survey 

■ 

CS| . A. Rawat^'^, Ajit K. Kembhavi^, F. Hammer^ H. Flores\ and S. Barway^ 

Oh" 
< 

I Context. Luminous Compact Galaxies (LCGs) (Mb < -20, R1/2 < 4.5 kpc and £^Wo(OII) >15A) constitute one of the most rapidly evolving 
O galaxy populations over the last ~8 Gyr history of the universe. Due to their inherently compact sizes, any detailed quantitative analysis of their 
^ morphologies has proved to be difficult in the past. Hence, the morphologies and thereby the local counterparts of these enigmatic sources have 
^ been hotly debated. 

1 I Aims. Our aim is to use the high angular resolution, deep, multiband HST/ACS imaging data, from the HST/ACS GOODS survey, to study the 

quantitative morphology of a complete sample of LCGs in the redshift range 0.5 < z < 1.2. 

Methods. We have derived structural parameters for a representative sample of 39 LCGs selected from the GOODS-S HST/ACS field, using 
^ full 2-dimensional surface brightness profile fitting of the galaxy images in each of the four filters available. ^435^^ - Zssolp color maps are 
. constructed for the sample to aid in the morphological classification. We then use the rest frame B band bulge flux fraction (B/T) to determine 
the morphological class of galaxies which are well fit by a bulge+disk two dimensional structure. Mergers were essentially identified visually 
by the presence of multiple maxima of comparable intensity in the rest frame B band images, aided by the color maps to distinguish them from 
HII regions. We also make use of the Spitzer 24yL/m source catalog of sources in the CDFS to derive the dust enshrouded star formation rates 
(SFR) for some of the sample LCGs 

Results. We derive the following morphological mix for our sample of intermediate redshift LCGs: 

Mergers: ~36%, Disk dominated: ~22%, SO: ~20%, Early types: ~7%, Irr/tadpole: ~15%. We establish that our sample LCGs are intermediate 
' mass objects with stellar mass ranging from 9.44 < Logio(M/MQ) < 10.96, with a median mass of Logio(M/MQ) = 10.32. We also derive SFR 
^ values ranging from a few to ~ 65 Mo/year as expected for this class of objects. We find that LCGs account for ~26% of the Mb < -20 galaxy 
kJ3I population in the redshift range 0.5 < z < 1.2. We estimate a factor ~11 fall in the comoving number density of blue LCGs from redshifts 
rN 0.5 < z < 1.2 to the current epoch, even though this number is subject to large uncertainities given the small sample size at zero redshift 
available from the literature. 

Conclusions. The strong redshift evolution exhibited by LCGs, and the fact that a significant fraction of LCGs are in merging systems, seem to 
indicate that LCGs might be an important phase in the hierarchical evolution of galaxies. We envisage that some of the LCGs that are classified 
as merging systems, might go on to rebuild their disks and evolve into disk galaxies in the local universe. 

Key words. Galaxies: starburst - Galaxies: fundamental parameters - Galaxies: formation - Galaxies: evolution 
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1. Introduction 

Understanding the true nature of Luminous Compact Galaxies 
(LCGs) is one of the outstanding problems in observational ex- 
tragalactic astronomy. LCGs constitute one of the most rapidly 
evolving galaxy populations at intermediate redshifts, with 
up to a factor ~10 fall in the comoving number density from 
a redshift of ~1 to the present epoch (Phillips et al. 119971 
Werk et al. 12004 ). They were first reported by Koo and Kron 
(119881) as possible blue compact galaxies due to the presence 
of strong narrow emission lines in the spectra of point-like 
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sources, which were otherwise believed to be QSO candidates. 
Hubble Space Telescope (HST) follow up observations (Koo 
et al. 11994b of these compact narrow emission line galaxies 
(CNELGs) gave typical half light radii of -0.65 arcsec. The ob- 
servations also provided evidence for exponential light profiles 
for these objects (as opposed to r^^^ profiles). Using high reso- 
lution spectroscopy from HIRES on the Keck, Koo et al. (1995} 
reported that these CNELGs have roughly Gaussian profile, 
narrow emission lines with velocity widths ranging between 
(T = 28 - 157 km s"^ They suggested that these CNELGs might 
be the progenitors of local spheroidal galaxies by fading of upto 
4-5 mags. 
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Guzman et al. ([T9971) and Phillips et al. ([T9971) identified 
a large population of such compact galaxies, defined as those 
with half-light radii ri/2 < O.Sarcsec and high surface bright- 
ness within half light radius (yU/8i4 < 22.2 mag arcsec"^) in the 
flanking fields of the Hubble Deep Field (HDF). Using high 
resolution Keck spectroscopy, they determined the emission 
line profiles to be roughly Gaussian with velocity widths rang- 
ing from cr ~ 35kms"^ to 150 km s"^ However, they stated 
that most compact galaxies at moderate redshifts yielded little 
morphological information even in HST/WFPC2 images. They 
reported about 2/3 of the galaxies in their sample to be consis- 
tent with being young star forming HII galaxies. 

Guzman et al. ([T998b found in HST/WFPC images of five 
CNELGs the presence of blue high surface brightness knots 
surrounded by a diff'use "exponential like" component. The 
knots were identified as the location of the current star forma- 
tion due to their bluer colors, whereas the diff'use exponential 
like component was interpreted to be an older underlying pop- 
ulation. 

More recently Hammer et al. (1200 lb have used VLT spec- 
tra of LCGs to argue in favour of their being the progenitors 
of present day spiral bulges, with tidally pulled in gas from in- 
teracting systems fueling high rates of star formation, resulting 
in the completion of the bulge formation and beginning of the 
disk formation. 

Hammer et al. (120051) have used multiwavelength observa- 
tions of 195, z > 0.4 intermediate mass galaxies to propose a 
scenario in which a majority of intermediate mass spirals have 
experienced their last major merger event (leading to the dis- 
ruption of the disk) within the last 8 Gyr. The merger scenario 
seems obvious since the merger rate was about 10 times more at 
z~l than it is today (Le Fevre et al. 12000 ). This merger phase is 
followed by enhanced star formation due to the merging, which 
is termed as the LCG phase, and a subsequent inside-out re- 
building of the disk leading finally to the formation of grand 
design spirals that we see today. Hence the LCGs are being 
proposed as a phase in the evolutionary history of galaxies in 
the hierarchical scenario of galaxy formation. 

Noeske et al. ('2006') have analyzed the morphologies of 26 
Luminous Compact Blue Galaxies (LCBGs) in the HST/ACS 
Ultra Deep Field atz~0.2 - 1.3. They concluded that the ma- 
jority of high z LCBGs are small galaxies that will evolve into 
small disk galaxies and low mass spheroidals in the local uni- 
verse. However, the sample selection of Noeske et al. includes 
galaxies as faint as Mb < -18.5 and, as per their own opinion, 
is not fully comparable with the sample selection criterion of 
Hammer et al. (120011) (which might favor progenitors of larger 
local galaxies), making any comparison of their results rather 
uncertain. 

In this paper we have made use of the high resolution im- 
ages provided by the Advanced Camera for Surveys (ACS) 
on board the HST to determine the quantitative morphology 
of intermediate redshift LCGs using the public data from the 
HST/ACS GOODS survey (GiavaHsco et al. 2004) . The unique 
advantage provided by the GOODS survey is that ACS imag- 
ing (as opposed to WFPC2 imaging) with it's improved driz- 
zled pixel scale of 0.03 arcsec/pixel is a factor of ~2 better than 
the typical drizzled pixel scale of WFPC2. This better sam- 



pling becomes specially important for compact sources such as 
LCGs with typical ri/2 < 0.5arcsec, for which the whole ex- 
tended emission (eg. disk) may be sampled by just a few pixels 
of the WFPC2. This in turn allows us to classify galaxy mor- 
phology on the basis of their luminosity profile, even for com- 
pact sources. In addition to this, the large area coverage of the 
GOODS survey provides a large enough sample of LCGs to 
derive robust statistics of their properties. 

This paper is organized as follows: In Section 2, we briefly 
describe our sample and in Section 3 present the methodol- 
ogy of our morphological classification and basic results. In 
Section 4 we describe the use of data at other wavelengths 
available in the Chandra Deep Field South (CDFS). In Section 
5 we provide individual description for each of the 39 LCGs. 
In Section 6, we calculate the evolution in number density of 
LCGs from intermediate redshifts to the current epoch. Our re- 
sults are discussed in Section 7 and our conclusions summa- 
rized in Section 8. Throughout this paper we adopt a cosmol- 
ogy with Hq =70 km sec"^ Mpc"\ D.m = 0.3 and Qa = 0.7 un- 
less specially stated. 



2. The Data 

We have used in our work: 

1. The publicly available version vl.O of the reduced, cali- 
brated, stacked and mosaiced images of the CDFS acquired 
with HST and ACS as part of the Great Observatories 
Origins Deep Survey, GOODS (GiavaHsco et al. l2004 D. 

2. The publicly available SExtractor (Berlin and Arnouts 
119961) based version rl . 1 of the ACS multi-band source cat- 
alogs as released by the GOODS team. 

3. Spectroscopic redshifts taken from the publicly available 
redshift catalog of the Vimos VLT Deep Survey VVDsQ 
(Le Fevre et al. 2004) . and 

4. The near-IR J & Ks band imaging data from the 
ESO G OODS/EIS Release Version 1.0 (Vandame et al, 
in preparation) which was obtained as part of the GOODS 
using the ISAAC instrument mounted at the Antu Unit 
Telescope of the VLT at ESO's Paranal Observatory, Chile. 
This data release includes 21 fully reduced VLT/ISAAC 
fields in J and Ks bands, covering 131arcmin^ of the 
GOODS/CDFS region. 

The near-IR data, along with the four band HST imag- 
ing data, combined with the spectroscopic redshifts men- 
tioned above, were used to derive the rest frame B band 
absolute magnitudes Mb of the galaxies. A detailed discus- 
sion of how the K-corrections were estimated and and how 
we finally arrived at the absolute magnitudes is provided in 
Appendix A of this paper. 



^ We have not applied any explicit quality cut on the VYDS red- 
shifts. However due to the requirement of the presence of Oil emis- 
sion line for an LCG, we have a sample of objects with very robust 
redshift measurements, with ~2/3 of the objects having redshift with 
better than 95% confidence, and all but 2 objects having redshifts with 
better than 75% confidence. 
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2.1. Sample selection 

In order to select a sample of Luminous Compact Galaxies, we 
have used the three criteria defined by Hammer et al. (l2QQlb 

First, we apply a luminosity criterion, selecting only galax- 
ies with Mb(AB) < -20. In addition to justifying the class 
as Luminous, this also ensures that there is sufficiently large 
signal-to-noise ratio (SNR) in the observed images to allow a 
detailed bulge-disk decomposition of galaxies even in the high- 
est redshift bin. 

Second, a compactness criterion was applied using a com- 
pactness parameter 



Redshift and Mg distribution for our sample 



Sm^ = m^(4.5 kpc) - m,(13.5 kpc) < 0.75 



(1) 



This ensures that at least half the flux is contained within an 
aperture of radius 4.5 kpc. The size of the apertures used by us 
are taken to be equivalent to those of Hammer et al., after cor- 
recting for the different cosmology. This correction is applied at 
the median redshift of our galaxy sample. One important differ- 
ence is that we have used elliptical instead of circular apertures 
for the galaxies, using the THETA JMAGE and ELLIPTICITY 
parameters of SExtractor to generate unique elliptical apertures 
for each object. Hence, we are using the most flexible aper- 
tures possible which allow for all possible isophotal shapes, 
and do not make any assumptions regarding the shapes of the 
isophotes to be expected. The radius of the aperture mentioned 
above is therefore understood to mean the length of the semi- 
major axis of the ellipse. 

Third, we include only those galaxies with known [Oil] 
A3121 emission, with rest frame equivalent width [EWqIOII] > 
15A]. This requirement of the presence of [Oil] A3121 line in 
the observed spectrum puts a lower limit on the redshift of z 
~0.5 for a galaxy in our sample, as the minimum observed 
frame wavelength in the VVDS is -5600 A. We have also put 
an upper redshift cutoff of z=L2 as there are few sources be- 
yond this redshift in our catalog and incompleteness becomes 
a serious issue at higher redshifts. 

Applying these three selection criteria to the GOODS- 
S HST/ACS field yields a sample of 39 Luminous Compact 
Galaxies. This accounts for -26% of the Mb < -20 galaxy 
population in the redshift range 0.5 < z < 1.2. 

Figure [T] shows the distribution of the redshift and the ab- 
solute B band magnitude of the LCG sample. Table [T] lists the 
three parameters that have been used in the selection criterion 
for each galaxy, along with their redshift. The ID given in col- 
umn (1) is used by us throughout this paper, and the GOODS 
ID given in column (2) is as listed in the rL 1 of the ACS multi- 
band GOODS source catalogs. 
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Fig. 1. The histogram of the redshift and the absolute B band 
magnitude distribution for our LCG sample. 



the magnitude- size plane and plotted our 39 LCGs on the same 
graph. This is shown in Fig. O Our galaxies occupy a region 
on the bottom left of the graph as expected by high surface 
brightness objects. In this region of the magnitude- size plane, 
the GOODS survey is at least 90% completqj. This proves that 
our sample is not biased by incompleteness issues. 
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2.2. Completeness of the LCG sample 

In a flux limited survey like GOODS, one has to be careful to 
avoid bias caused by incompleteness of the sample, especially 
when one is working with faint objects. Since we are working 
with intrinsically luminous and compact sources, they are ex- 
pected to be high surface brightness objects free of such incom- 
pleteness biases. In order to demonstrate this, we used the com- 
pleteness curves published by the GOODS team as depicted on 



Fig. 2. F850LP magnitudes and half light radii (R1/2) of our 
39 LCGs shown in comparison with the completeness curves 
of the GOODS survey as depicted on the magnitude- size 
plane. The percentage completeness is marked on the individ- 
ual curves. The completeness curves are plotted in intervals of 
10% and were kindly provided by Harry Ferguson. 

^ In order to see how the completeness curves were estimated, refer 
to Giavalisco et al. (120041) 
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Table 1. The LCG sample 



Our 


GOODS 




Mb 


6m^ 


EWoiOin 


ID 


ID 


Redshift 


(Mag.) 


(Mag.) 


(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


904260 


J03321 1.97-275033.9 


0.983 


-20.73 


0.72 


35.39 


904604 


J033212.62-275105.6 


0.990 


-20.38 


0.51 


29.38 


904680 


J033212.78-275012.8 


0.964 


-20.86 


0.64 


70.19 


905632 


J033214.68-274337.1 


0.976 


-20.04 


0.48 


107.80 


905983 


J033215. 36-274506.9 


0.860 


-20.74 


0.49 


49.41 


906961 


J033217. 14-274303.3 


0.566 


-20.93 


0.29 


16.01 


907047 


J033217.28-275416.6 


1.112 


-21.23 


0.64 


23.73 


907361 


J033217.77-274714.9 


0.731 


-20.17 


0.60 


19.30 


907794 


J033218.52-275508.3 


1.144 


-21.44 


0.58 


46.77 


908243 


J033219.32-274514.0 


0.726 


-20.35 


0.73 


43.56 


909015 


J033220.70-275501.7 


1.039 


-20.25 


0.46 


49.15 


909093 


J033220.86-275405.1 


0.968 


-20.32 


0.62 


95.48 


909429 


J033221.43-274901.8 


0.737 


-20.81 


0.69 


33.67 


910413 


J033223.03-275452.3 


0.655 


-21.30 


0.54 


18.38 


911747 


J033225.20-275 100.1 


0.840 


-20.73 


0.47 


51.07 


911780 


J033225.26-274524.0 


0.664 


-20.88 


0.60 


20.71 


911843 


J033225. 35-274502.8 


0.973 


-20.49 


0.34 


67.75 


912744 


J033226.61-275001.9 


0.690 


-20.42 


0.60 


59.31 


913482 


J033227.62-274144.9 


0.664 


-20.69 


0.28 


45.17 


914038 


J033228.48-274826.6 


0.667 


-20.17 


0.43 


30.06 


915400 


J033230.43-275 140.4 


0.764 


-20.08 


0.67 


99.46 


916137 


J033231.50-275241.4 


0.980 


-20.14 


0.58 


50.62 


916446 


J033231.96-275553.4 


0.839 


-20.03 


0.60 


51.10 


916866 


J033232.61-275316.7 


0.987 


-20.06 


0.64 


29.44 


918147 


J033234. 61-275324.5 


1.099 


-21.52 


0.65 


71.64 


919573 


J033236.72-274406.4 


0.665 


-20.05 


0.69 


37.62 


919595 


J033236.74-275206.9 


0.785 


-20.40 


0.49 


49.14 


920435 


J033238.08-275248.7 


1.034 


-20.15 


0.43 


86.23 


921406 


J033239.65-275226.2 


1.095 


-21.12 


0.72 


23.93 


922675 


J033241.88-274853.9 


0.666 


-20.31 


0.54 


27.08 


922733 


J033242.02-275226.0 


0.650 


-20.04 


0.40 


55.70 


922761 


J033242.09-275 109.4 


0.961 


-20.10 


0.62 


35.14 


923085 


J033242.73-275 159.2 


1.122 


-21.12 


0.50 


47.24 


923926 


J033244.45-274940.2 


1.012 


-20.86 


0.38 


93.41 


924881 


J033246.59-274516.3 


0.839 


-20.05 


0.52 


68.14 


926109 


J033249.53-274630.0 


0.522 


-20.13 


0.47 


29.64 


926217 


J033249.87-275 129.0 


0.767 


-20.25 


0.64 


28.37 


907305 


J033217.68-274208.9 


1.185 


-20.76 


0.45 


68.82 


914895 


J033229.7 1-274507.2 


0.736 


-20.16 


0.56 


43.32 



6m^ is the aperture magnitude difference given in Equ0in section 2.1 



3. Morphologies of LCGs 

Using the HST/ACS GOODS images, quantitative two dimen- 
sional bulge-disk decomposition has been carried out to derive 
the structural parameters of the galaxies which are then used 
to quantify their morphologies. Information derived from the 
color maps is also used in the morphological classification. 



3. 1 . Structural parameters 



We used the software galfit (Peng et al."2002) to carry out two- 
dimensional modeling for our sample of 39 LCGs in each of 
the four HST/ACS filters to obtain their structural parameters. 
Galfit models the galaxy image as a linear combination of a 
bulge and a disk, using a well established analytic model for 
each of the two components. An important caveat in this ap- 
proach is the underlying assumption that a galaxy can be repre- 
sented as a linear combination of simple, smooth analytic func- 
tions. Real galaxies are known to be much more complex, with 
the presence of spiral arms, bars and central point sources etc. 
Despite these reservations, galaxy fitting algorithms have been 
demonstrated to be successful in the past for the purpose of 



quantitative morphological classification of galaxies (Peng et 
al. I2Q02I Simard et al. l2002l) . 

The intensity profile of the bulge is modeled with the Sersic 
law (Sersic HSMl), 



X(r) = S^exp 



-2.303/?, 



1 



(2) 



with re being the half light radius, 
re and n being the sersic index. 

The disk is modeled as an 
(Freeman 119701), 



the surface brightness at 
exponential function 



i:(r) = Soexpf-— 1 



(3) 



with So being the central surface brightness and the e-folding 
scale length. 

The model galaxy is then constructed as a linear combi- 
nation of these two components plus other free parameters 
such as the sky background value (bkg)x,y in pixel position 
(x,y), which is assumed to be constant across the galaxy im- 
age, the ellipticities etuige and Edisk of the bulge and the disk 
isophotes respectively and their position angles Otuige and Odisk- 
The SExtractor based photometric catalog of our sources pro- 
vided us with valuable starting values for parameters such as 
sky background, half light radii of the objects and their elliptic- 
ity and position angles etc. This analytic model of the galaxy is 
then convolved with the PSF of the observation and compared 
with the observed galaxy image. The values of the free param- 
eters are determined iteratively by minimizing the reduced 
defined as usual. 



1 (galaxy^ ^, - mode^3,) 



(4) 



where A^dof is the number of degrees of freedom, galaxy x,y and 
modelx,y are the counts in the pixel (x,y) of the galaxy and the 
model image respectively, cTx^ is the noise in the pixel (x,y) 
and nx,ny are the number of pixels in the galaxy image in the x 
and y direction respectively. 

Special attention is paid in arriving at a reliable noise model 
crx,y for use in Eq. ffl In case of the HST/ACS GOODS fields, 
the noise characteristics of the science images are stored in 
the form of weight maps, Wx,y which are images produced as 
part of the data reduction process called drizzling (Fruchter 
& Hook 1997120021) , which give a measure of the background 
plus instrumental noise per pixel in the science data. Since the 
GOODS images are mosaics consisting of a diff'erent number 
of overlapping pointings at diff'erent areas, the eff'ective expo- 
sure time varies across the field leading to variable depth across 
the field. The weight maps take into account this variation and 
should provide a faithful estimation of the inverse variance 
(l/crly) per pixel. However, the interpolations introduced by 
drizzling result in correlations between pixels in the drizzled 
science images. The weight maps are normalized to show the 
expected noise per pixel that the images would have in the ab- 
sence of these correlations. Therefore, the apparent cTxj that 
one measures in the science image is smaller than that given by 
the (1/Wjcj)^^^, because the apparent cTx^y is suppressed by the 
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effects of these correlations. This necessitates the use of a cor- 
rective factor when deriving the cr^^ y from the weight maps Wx,y 
such that: 



Fx 



1 



(5) 



where F depends upon the drizzling parameters SCALE (s) 
and PIXFRAC (p) (see Fruchter & Hook 1 199712003) . For the 
GOODS dataset, s=0.6 for all the four filters, whereas p=0.8 
for F435W, F606W and F775 filters, and p=0.7 for F850LP 
filter (Anton Koekemoer, private communication). Using the 
prescription given in Casertano et al. (12000b , we derive a value 
of F=0.61 for the F850LP filter and F=0.56 for the other three 
filters. These derived values of F, in conjunction with Eq. \5\ 
and the supplied weight maps provided the noise model for our 
sample of galaxies. A more detailed discussion on weight map 
conventions and noise correlation in drizzling can be found in 
Casertano et al. (2000). 

We have determined the sky background using SExtractor 
and held it fixed at this value during the minimization. All other 
parameters like r^, r^, n, position angle etc. were allowed to 
be free in the fitting process, so that galfit can explore the full 
range of parameter space to find the global minimum and is not 
pulled into a local minimum due to constraints put in by hand. 

Once the best fit galaxy model has been obtained, we cal- 
culate the bulge fraction B/T, defined as: 



B 
f 



flUXjyulgg 



fluXbulge + flu^disk 



(6) 



where the flux^uige & flu^disk are calculated by integrating the 
bulge and the disk profiles respectively over all values of the 
semi-major axis using the best fit parameters obtained by gal- 
fit. This bulge fraction B/T is found to be broadly correlated 
with the traditional Hubble type of a galaxy in the sense that an 
early type galaxy has a high B/T ratio and a late type spiral has 
a low value. It is this ratio that we use for characterizing the 
morphology of our galaxy sample. 

In the above scheme, an accurate estimation of the PSF is 
crucial to the robust determination of the derived structural pa- 
rameters of a galaxy. In our work, the star closest to a given 
galaxy is used to determine the PSF for the purpose of bulge- 
disk decomposition. Stars are the most pristine characterization 
of the PSF in any given observation, as they have gone through 
the same optics as well as the reduction pipeline as the pro- 
gram galaxies. Only bright (non saturated) and well behaved 
stars, free of neighbours and other contamination were used as 
candidate PSF stars in our work. A master list of such stars 
was prepared by visually examining all candidate stars in the 
GOODS field and rejecting any unsuitable star based on the 
basic criterion outlined above. This yielded ~60 stars scattered 
all across the CDFS which could be used to determine the PSF, 
from among -400 candidate stars. 

We looked for the variation of the PSF half-light radius with 
RA/DEC for the -60 PSF stars that we had shortHsted. A "sur- 
face plot" of the PSF half light radius variation is shown in 
Fig. [3l The PSF is seen to be extremely constant across the 
GOODS-S field of view. Barring a couple of "peaks" in the 



plot (which are explained on the basis of the red halo eff'ect be- 
low), the 7? 1/2 of the PSF stars is stable to within 5%. This is 
because of the extremely stable detector/optical system off'ered 
by HST/ACS which yields almost identical PSF even for two 
images taken at diff'erent epochs or for two diff'erent exposures 
(as in GOODS). In addition, the lack of an atmosphere takes 
away a major source of variability in the PSF from one epoch 
to another seen from the ground. 

PSF variation across the CDFS in the z-band 




53.09 

RA{degrees) 



Fig. 3. The variation of half-light radius(in units of ACS pixels) 
of stars with RA and DEC across the GOODS-S field in the 
F850LP filter. This plot uses ~60 stars spread across the field. 

One important point to consider while interpreting Fig. [3] 
is that whatever variation of PSF is seen in this figure is due 
to the combined eff'ects of both spatial variation of the PSF, 
plus any other factors such as the red halo (Sirianni et al. 120051 
GilHland and RiessTOO?) eff'ect in the F850LP band ACS PSF 
which depends on the color of the star. In fact, we find that all 
of the "peaks" that are seen in Fig. [3] belong to stars that are 
extremely red, and therefore the red halo might be responsible 
for their large 7? 1/2 rather than any spatial variatiorH. 

Extensive tests were carried out to rule out the possibility 
of variations in the determination of the structural parameters 
of galaxies due to the use of diff'erent stars as the PSF. For a 
sample of -450 galaxies in the CDFS, we used galfit in batch 
mode to derive the structural parameters and B/T ratio twice 
for each galaxy. In the first run, we used the star nearest to the 
galaxy being modeled as the PSF. In the second run, we used 
the second nearest star to the galaxy as the PSF, everything 
else remaining the same as in the first rur0. In Fig.|4]we show 
the histogram of the diff'erence between the B/T ratio obtained 
in the two cases. As is evident from the figure, the estimation 
of B/T ratio is quite robust irrespective of which star is used 
as the PSF. About 70% of the sources have their B/T value 



^ For a more detailed discussion of the red halo effect on the 
F850LP PSF stars used by us and how we dealt with it, see Appendix 
B 

^ The median distance from the sample galaxy to the nearest PSF 
star is 48 arcsec while the median distance to the second nearest PSF 
star is 80 arcsec. 
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stable to within 10%. The sharp peak around zero means that 
the B/T ratios obtained in the two run agree quite well, and the 
symmetrical shape of the peak at zero points towards random 
errors as the cause for any dispersion in the B/T values obtained 
in the two runs rather than any systematic effect. 

We conclude that the final value of B/T is stable to within 
10% irrespective of which star is used as the PSF. Keeping this 
in mind, the star closest to a given galaxy is used as the PSF 
for the purpose of bulge-disk decomposition of our sample of 
LCGs. The important thing to notice is that in Fig.lH the diff'er- 
ences in B/T ratio that we obtain on using two diff'erent PSFs is 
due to the cumulative eff'ects of PSF variation due to all possi- 
ble reasons such as spatial variation of the PSF, red halo eff'ect 
etc. Since the cumulative error in the derived B/T ratio is shown 
to be within ~10%, our results are stable to within this error ir- 
respective of the fact as to which particular factor contributes 
to the variation in the PSF. This should convince the reader that 
our final results are insensitive to the choice of PSF star. 

Figure \5\ shows the result of our quantitative bulge-disk 
decomposition for one of the LCG sample galaxies at a red- 
shift of 0.976. The upper left frame shows a 30kpc x 30kpc 
F850LP galaxy image with the galaxy name, redshift and 
morphological classification indicated. The upper right frame 
shows the best fit model image having ;^^^^= 1.09 3. Notice that 
the galaxy in the neighbourhood of the primary target (north- 
east) is also being fitted simultaneously. In this work, we have 
taken the approach of fitting a single component sersic func- 
tion to all neighbouring objects within ~5.0arcsec of the pro- 
gram galaxy. This approach, though more computation inten- 
sive, gives vastly improved results compared to the traditional 
approach of masking the neighbouring objects in the frame 
(which does nothing to account for the flux contributed to the 
program object by the extended wings of the neighbours). This 
efl'ect can be particularly severe in case of bright neighbours 
and would lead to an inaccurate estimation of the structural pa- 



Fig. 5. From top-left to bottom-right 30 x 30 kpc F850LP im- 
age, Galfit Model image. Residual image and the B-z color map 
of the LCG 905632. North is up. East is to the left. The galaxy 
name, redshift and morphological classification are indicated 
in the F850LP galaxy image. The horizontal bar at the bottom 
right of the galaxy image indicates 1 arcsec. The intensity scal- 
ing is the same (squareroot) in each of the three grayscale im- 
ages. 

rameters of the program galaxy. At the bottom left, the residual 
image is shown, which is the diff'erence between the galaxy im- 
age and the best fit model. The residual image is devoid of any 
structures and consists solely of Poisson noise, which is an ad- 
ditional indication of the good quality of the fitting. Finally, the 
bottom right frame shows the B-z color map of the galaxy in 
AB magnitude system; a blue core is seen. This is a disk dom- 
inated galaxy with a B/T ratio of 0.33 and a disk scale length 
rd = 1.16 ±0.02 kpc. 

Figure[6] shows a comparison between the azimuthally aver- 
aged radial surface brightness profile of the same LCG 905632 
(points with error bars), and the surface brightness profile of the 
model image produced by galfit (red smooth line), produced 
using the IRAlU task ELLIPSE. The two profiles agree within 
errorbars, which serves as an additional test of the accuracy of 
the fit. 

3.2. Morphological classification 

We have carried out bulge-disk decomposition on all the 39 
LCGs in our sample. The galaxies that are well fit by a 
bulge-hdisk two-dimensional structure, are classified into three 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc. under cooperative agreement with the 
National Science Foundation. 
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Fig. 6. A comparison of the azimuthally averaged F850LP band 
radial surface brightness profile of the galaxy 905632(points 
with error bars), as compared to the surface brightness profile 
of the model image produced by galfit(red smooth line) pro- 
duced using the IRAF task ELLIPSE. 



'Hubble types' mainly based on the rest frame B band BjT 
ratio: early types E (0.8 < BIT < 1), intermediate types SO 
(0.4 < BIT <0.8) and disk dominated S (0.0 < BIT <0.4). 
This segregation of galaxies into diff'erent Hubble types based 
on their B/T ratio is made using the criterion described by 
Zheng et al. (120041 ^05 ), even though we have grouped our 
galaxies more broadly in B/T as compared to Zheng et al. Since 
our sample of LCGs range in redshift from 0.5 to L2, the rest 
frame B band falls in diff'erent observed filter bandpasses, de- 
pending upon the redshift of the object. We made use of the 
multiband nature of the HST/ACS GOODS dataset to apply a 
morphological k-correction. What this means is that we use the 
structural parameters from the observed frame band that corre- 
sponds most closely to the rest frame B band, given the redshift 
of the object. For this, the redshift range of the galaxies and the 
observed frame filter used for obtaining the structural parame- 
ters are paired as: 

0.5 < z < 0.6 (F606W); 0.6 < z < 0.85 (F775W); 0.85 <z< 
1.2 (F850LP). Since the minimum redshift of a galaxy in our 
sample is -0.5, the observed frame F435W fits were not used in 
our classification scheme. The structural parameters obtained 
by galfit for galaxies which are well fit by bulge-hdisk two- 
dimensional structure in F850LP, F775W and F606W bands 
are given in Table[2l[3]and|4]respectively. We find that galaxies 
that have redshifts at the border between any two of the three 
redshift bins described above have very comparable structural 
parameters in both the filters. We are unable to obtain accept- 
able fits for some galaxies in the shorter wavelength filters even 
though they have reasonable fits in the F850LP band. This is 



primarily due to the fact that at intermediate redshifts, the bluer 
filters are probing the UV region of galaxies which have non- 
uniform and clumpy light distribution, making it difficult to ob- 
tain a meaningful fit. 

In addition to the three 'Hubble types' mentioned above, 
an additional type, tadpole/irregular (Irr) is included in our 
classification scheme to describe the objects without a clear 
bulge-hdisk structure. A single nucleus off'set from the "cen- 
ter" of the image (i.e. having a one sided elongated light dis- 
tribution) prompted us to tag some objects as "tadpole". We 
designated a galaxy as irregular if it had a patchy, blue light 
distribution with no clear nucleus at all. 

A quality factor running from 1 (secure fit) to 4 (fit failed) 
is used to give the confidence of the galfit fit. The quality flags 
were decided using the following guidelines: 
Q=l: Perfect fit Wiihxled - ^-^ bulge mag error < 0.1 and 
disk mag error < 0.1 and fractional error in bulge effective ra- 
dius Tg, disk scale length and Sersic index n < 10% 
Q=2: One or two of the above conditions are violated or the 
presence of an AGN is indicated due to CXO X-ray detection 
which might complicate the bulge-disk fitting. 
Q=3: More than two of the above conditions are violated or a 
fit was overruled by visual inspection to classify the galaxy as 
Irr/Tadpole. 

Q=4: The fit fails completely due to the presence of multi- 
ple nuclei or patchy light distribution and fits the bright nuclei 
rather than the underlying component, or if the galaxy is clas- 
sified as a merger by visual inspection. 

In some cases, the above quality ffag had to be overruled 
where one of the components is very weak (eg. 905983), re- 
sulting in large uncertainties in the determination of parameters 
of the weaker component (bulge in case of 905983). This is in- 
dicative of the fact that there are not enough counts in that par- 
ticular component leading to large errorbars on the determined 
parameters. In such a case, the B/T ratio has large fractional 
errors, but it still provides a robust discrimination between a 
bulge dominated and a disk dominated system, and it is not 
proper to tag it as a poor fit. 

We have used the method of elimination to classify a galaxy 
as a merger. If a galaxy has a sufl&ciently smooth brightness 
profile allowing for a good fit to be obtained by galfit, it is clas- 
sified as E/SO/S depending on it's B/T ratio. If however galfit 
fails to provide a good fit (as per the criteria given above), it 
is supposed that the brightness distribution of the galaxy is too 
distorted to be put in either of the above three morphological 
classes. It is at this point that we invoke the use of colormaps 
and/or visual examination of ACS images to discriminate be- 
tween irregulars and mergers. Mergers were essentially identi- 
fied visually by the presence of multiple maxima of comparable 
intensity in the rest frame B band images. This was aided by the 
use of color maps to distinguish multiple nuclei (generally red- 
der) from HII regions (generally bluer). A label is employed to 
identify possible merging systems (M1/M2); see the footnotes 
of Table [5] for details. 

Even though in the first run, galfit was used in automated 
batch mode, the results were examined visually for obvious 
deficiencies such as badly fitted neighbours, or a badly fitted 
program galaxy due to patchy light distribution or presence of 



8 



A. Rawat et al.: Morphologies of Luminous Compact Galaxies 



bright HII regions in the galaxy. In many cases where bad fits 
were obtained, a decision had to be made, for example as to 
whether the patchy light distribution is due to the presence of 
HII regions, or due to the presence of double nuclei. Galfit, by 
design, cannot provide the necessary information to identify a 
merger. Hence it is not sensitive to subtle diff'erences between 
multiple nuclei and HII regions (for example) and visual ex- 
amination is needed in such cases. In order to avoid individual 
biases, in the classification of such objects, visual inspection of 
the galaxy's appearance in HST/ACS imaging and color distri- 
bution were performed by two of the authors (A.R. and S.B.) 
independently. A final classification was made after discussing 
any inconsistent cases. Table [5] lists the results of our morpho- 
logical classification, including galaxy type (Col. 9), quality 
factor (Col. 10) and merging type (Col. 12), and the individual 
description for each object is given in Section 5. 

3.3. The morphological mix of LCGs 

We have derived the following morphological mix for our sam- 
ple of intermediate redshift LCGs: 

Mergers: -36%, Disk dominated: -22%, SO: -20%, Early 
types: -7%, Irr/tadpole: -15%. 

The detailed results of our morphological classification are 
listed in Table [21 and individual comments for each galaxy are 
given in Section 5. The significance of the relative fractions of 
morphological types that we obtained is discussed in Section 7. 

3.4. Color maps 

B-z color maps were produced for our entire sample of 39 
LCGs. We directly used the F435W and F850LP vl.O images, 
which are astrometrically aligned with each other to an accu- 
racy of better than a fraction of a pixel. 30 kpc x 30 kpc cutouts 
of these color maps and the associated F850LP grayscale im- 
ages for each of the 39 LCGs are shown in Fig.[8l 

We also made use of the integrated observed colors of these 
galaxies by comparing them with the color-redshift evolution 
curves as predicted by Bruzual-Charlot stellar population syn- 
thesis code GALAXEV 2003 (Bruzual and Chariot 1993 2003) 
for diff'erent galaxy types. We have shown in Fig. [7] three dif- 
ferent evolution models, all for galaxies with Solar metallicity 
formed at redshift=5.0: an instantaneous burst corresponding to 
early types, an exponentially decaying SFR with an e-folding 
time scale of 1 Gyr corresponding to SO, and another with a 
time scale of 7 Gyr corresponding to late types. The e-folding 
timescales for diff'erent Hubble types are the same as those used 
by Zheng et al. ( 2004). The redshift evolution of their B-z ob- 
served color is plotted along with the known B-z color and red- 
shift of the 39 LCGs segregated by their morphological class. A 
young starburst has also been shown for comparison, assuming 
a power law spectrum Ly ocy~^. 

In Fig. [T] many LCGs lie close to the curve for late type 
galaxies. It is also seen that most galaxies classified as mergers 
are close to the late type curve (i.e. they are bluer), as opposed 
to E/SO classified galaxies which are much redder (even though 
there is a lot of scatter). Hence, our classification acts as a crude 
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Fig. 7. B-z observed color vs redshift evolution for three. Solar 
metallicity, formed at redshift=5.0 spectrophotometric galaxy 
models assuming instantaneous burst (early type), exponen- 
tially decaying SFR with an e-folding time scale of 1 Gyr (SO) 
and 7 Gyr (late type). A young starburst has also been plotted 
for comparison, assuming a power law spectra Ly oc y"^ Our 
sample LCGs are plotted using diff'erent symbols for diff'erent 
morphological class, with Early Types (E) as red open circles, 
SO as black diagonal crosses. Disk dominated (S) as blue stars, 
Irr/Tadpole as blue open triangles and Mergers(Ml/M2) as blue 
filled squares. 



spectrophotometric galaxy classification scheme in which the 
galaxies designated by us as E/SO are generally redder and lie 
closer to the Early type or SO curve, than those designated as 
disk dominated/mergers/Irr. Every one of the galaxies classified 
as disk dominated (S) by us is redder than the late type model 
curve. This might be due to the eff'ect of dust reddening. Also, 
some of them can be early type spirals, which are expected to 
be redder than the late type curve. Most LCGs are seen to be 
redder than what is expected of a young starburst which is bluer 
than a typical late type galaxy (as shown). This might be due to 
the presence of an older underlying stellar population as estab- 
lished by Hammer et al. (120011) . Also, one must remember that 
our B-z observed color for individual LCGs are not corrected 
for dust extinction, where as the Bruzual and Chariot gener- 
ated curves are dust free models. The high rate of detection in 
mid-IR Spitzer (8/39 LCGs are detected in Spitzer, a -20% de- 
tection rate) indicates that these are dust enshrouded starbursts 
where the UV light from young 0/B stars is being reprocessed 
by the dust and is being emitted in mid-IR. This might help to 
explain the redder than expected colors for some of the LCGs 
in our sample. 
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4. Input from other wavelengths 

We have looked for the detection of our sample LCGs in other 
wavelength bands, to calculate their stellar masses, star forma- 
tion rates (SFR) & to look for possible AGN activity etc. Since 
CDFS is one of the most data rich regions on the sky, we have 
access to public data on objects in this region at various wave- 
lengths such as radio. X-ray and mid-infrared. The results of 
our analysis in these wavelength bands are given below. 

4.1. Radio emission from LCGs 

Koekemoer et al. (I2QQ6I) have prepared a comprehensive cat- 
alog of all radio sources in an area of 1.2 deg^ covering the 
CDFS to a limiting (Icr) sensitivity of ~14yuJy at 1.4 GHz us- 
ing the Australia Telescope Compact Array (ATCA). Afonso 
et al. r2006) have published optical and X-ray identifications 
of faint radio sources in the GOODS-CDFS ACS field. From 
this catalog we find that only 1 out of 39 LCGs of our sam- 
ple had any radio emission. The radio flux of this sole de- 
tected galaxy 906961 is 79yuJy (see individual description of 
this galaxy in Section 5). This puts an upper limit of ~14yuJy at 
1.4 GHz for the radio flux from the other 38 LCGs in our sam- 
ple. Assuming a power law radio spectrum Ly oc y"^, with a 
spectral index a = 0.5 (which is typically used as the dividing 
line between steep and flat spectrum radio sources), and using 
the highest redshift for our LCG sample (z=1.2), puts a strict 
upper limit of ~ 10^^ erg sec" ^ on the 1.4 GHz luminosity for 
the 38 non-detected LCGs. Since the VVDS spectrum does not 
have enough spectral coverage to discriminate between AGN 
and starbursts using line ratios, we decided to use X-ray and/or 
radio detection as a sign of AGN activity, so as to make sure 
that our LCG sample does not sufl'er from any serious AGN 
contamination. The limit that we derived for the 1.4 GHz lu- 
minosity rules out any radio loud AGN contamination of our 
sample. 

It must be mentioned here that the beam size of the ATCA 
radio survey was rather broad at 16.8 arcsecx6.95 arcsec. Since 
the typical size of our galaxies is of the order of ~ 1 arcsec, this 
restricts the usefulness of this radio catalog for our purpose. 

4.2. CXO X-ray detection 

Giacconi et al. (2002") have published a X-ray source cat- 
alog obtained from the -IMsec exposure of the Chandra 
Deep Field South (CDFS), using the Advanced CCD Imaging 
Spectrometer (ACIS) on the Chandra X-ray Observatory, as 
well as their optical counterparts. 

We used a very simple prescription of looking for potential 
X-ray counterparts for our sample of 39 galaxies by searching 
for X-ray sources within a 6 arcsec search radius of the program 
galaxy using the NASA Extragalactic Database (NED). Each 
potential X-ray counterpart was then checked visually using 
optical counterpart image cutouts published by Giacconi et al. 
(120021) . 

This simple method resulted in X-ray emission being con- 
firmed from four of the 39 LCGs. All four X-ray sources are 
confirmed to be AGN by Zheng W. et al. (l2QQ4b . The hard- 



ness ratio (HR) for these galaxies as given by Giacconi et al. 
(I2QQ2I) is listed in Table \5\ These are mainly provided for the 
reader who might want to distinguish between Type 1 and Type 
2 AGNs, although such a distinction was not necessary for our 
work (see Zheng W. et al. l2QQ4b . Interestingly, we find that all 
the 4 LCGs with X-ray emission are E/SO type galaxies. Two of 
these LCGs show clear blue point like nucleus in the colormaps 
(Fig. [8]), while in the other two cases, the AGN might be ob- 
scured. We also noticed that 3 of these X-ray detected LCGs 
are detected by Spitzer MIPS (see next section). This is not 
surprising given the known AGN contribution to mid-infrared 
surveys (Fadda et al.i2002j). 

4.3. Spitzer MIPS detection and SFR for LCGs 

The GOODS Spitzer Legacy program observations cover two 
fields on the sky. One of these fields (GOODS-N) coincides 
with the Hubble Deep Field North, while the other (GOODS- 
S) coincides with the Chandra Deep Field South. 

We made use of the public data made available as part 
of The Great Observatories Origins Deep Survey (GOODS), 
Spitzer Legacy Data Products, Third Data Release (DR3). In 
addition to other data products, this release consists of the 
"best-eff'ort" reductions of 24 micron data for the southern 
GOODS field taken with the Multiband Imaging Photometry 
for Spitzer (MIPS), and a 24 micron v0.91 source list of all 
sources brighter than 80yuJy. Dickinson et al. ("2006") describe 
the Spitzer observations, data reductions, and the data products. 
Chary et al. (l2QQ6b explain the preparation of the MIPS catalog 
in detail. 

We performed a simple positional cross correlation of our 
39 galaxy sample with the v0.91 MIPS source list with a toler- 
ance radius of 2.5 arcsec, which was chosen keeping in mind 
the ~ 5 arcsec PSF of the Spitzer MIPS observations. This 
yielded 24yum fluxes for 8 of the 39 LCGs. It was found that the 
maximum separation between the HST- Spitzer counterparts in 
case of a true match was of the order of -0.5 arcsec. Visual 
examination of the images was employed to ensure a one-to- 
one correspondence between the HST and the Spitzer sources. 
The other 3 1 LCGs for which no counterpart was found in the 
Spitzer MIPS catalog, have an upper limit of 80yuJy for their 
24yum fluxes. We calculated the infra-red luminosities Ljr for 
the 8 galaxies using the approximations given in Chary and 
Elbaz (i2001b . using the rest frame 12 and 15 jim luminosi- 
ties, which were in turn determined using the 24yum MIPS 
fluxes. We used the approximation of 15yum for galaxies with 
0.5 < z < 0.8 and that of 12yum for galaxies with 0.8 < z < 1.2 
to calculate the IR luminosity. The IR luminosities obtained 
using the two approximations agree to within 0.1 dex. We also 
calculated the upper limits for the infra-red luminosities Ljr for 
the other 3 1 LCGs with no counterparts in Spitzer/MIPS cata- 
log. These values for the Ljr are listed in Table [H 

Kennicutt (1998) has transformed the IR luminosity of 
young (age < lO^yr) starburst galaxies to a star formation 
rate (SFR). He gives an approximate estimate of the dust- 
enshrouded SFR p 

piM^yr-') = 1.71 X 10-i%^(Lo) (7) 
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Using this relation, we obtain SFRs varying from 20 - 
65Mojr"^ for our sample LCGs with 24yum MIPS fluxes, 
which is the range expected for galaxies of this class (eg. 
Hammer et al. 1200 ID . For the objects which have no detection 
in Spitzer/MIPS catalog, this method yielded an upper limit on 
their SFR. These values are listed in Tabled 

In order to be able to better constrain the SFR of those 
objects that are not detected in the Spitzer/MIPS catalog, we 
used their rest frame UV continuum flux to derive an estimate 
of the SFR. This method assumes that the rest frame UV flux 
that is observed for the objects is directly emitted by the 0/B 
type stars in a starburst and there is no dust extinction (see 
Kennicutt 1 19981 for details). Since galaxies are dusty, the UV 
flux undergoes varying amount of extinction, and hence the 
SFR yielded by this method is to be treated as a lower limit. 
Hence the SFR derived by this method are listed as lower lim- 
its in Table [3 

4.4. Stellar masses of our sample LCGs 

Stellar masses for our sample LCGs were calculated on the ba- 
sis of their absolute K band magnitudes, corrected for the pres- 
ence of red supergiant stars, using the empirical relationship 
of Bell et al. ( 2003i). A variation of the Salpeter Initial Mass 
Function (Salpeter 11955b , called diet Salpeter IMF was used 
as explained in Bell et al. and the values were calculated in 
a consistent manner as in Hammer et al. (2005 ). The derived 
values of stellar mass for our sample are listed in Table \5\ To 
derive the values of stellar mass using a classical Salpeter IMF, 
it is necessary to add 0.15 dex to the stellar mass values listed 
in Tabled Accuracy is better than 0.1-0.2 dex as it has been 
shown to correlate particularly well (for rotating disks) with 
log(Vm«x) at z=0.6 by Flores et al. (12006b . 

5. Individual description 

In Fig. [51 color map stamps (right) of the 39 LCGs are shown, 
along with F^SOLP imaging gray level stamps (left). The color 
bar in each color map stamp shows the B-z color range from 
to 4, except for objects ID 904604, 907361 and 909015 
where it ranges from to 6. A description of each target 
is presented below. The FITS files of all our 39 LCGs in 
each of the four filters can be downloaded from the follow- 
ing url http://ww w.iucaa.ernet.in/~ rawat/lcg -l .html[ They are 
far more informative than the postscript versions depicted here. 

904260 This is a possible ongoing merger. Spiral arms are dis- 
cernible in one of the constituent galaxies, specially in the B 
band image. The nucleus of the other member is extremely 
blue indicating star formation. The colormap also shows the 
spiral arms to be bluer than the rest of the galaxy. 

904604 This is a featureless edge on disk galaxy. It practically 
drops out of the blue filter and is hence quite a red object 
with a rest frame (U - B)o = 1.47 mag. 

904680 This object has two nuclei which becomes apparent in 
the shorter wavelength images as well as in the color map. 
Both the nuclei are bluer compared to the surrounding area. 
This is possibly an ongoing merger. 



905632 This is a relatively featureless disk galaxy. The col- 
ormap shows that the core is bluer compared to the outer 
edges of the galaxy. Also the Oil emission line is very 
strong at 107. 8 A. However, there is no detectable X- 
ray/radio emission from this galaxy as described in the text. 
Also there are no signs of nuclear activity, hence there are 
no indications for an AGN. 

905983 This looks like a face on spiral galaxy. The spiral arms 
are easily visible in the summed image as well as in the gal- 
fit residuals, especially in the shorter wavelength fits. Bright 
HII regions can be seen embedded inside the spiral arms in 
B and V band images. The colormap shows that the spiral 
arms are bluer than the rest of the galaxy. 

906961 This has the appearance of an early type galaxy with 
a distinct blue point like core as shown by the color map. 
A satisfactory bulge-disk decomposition could not be ob- 
tained due to the presence of the strong point source at the 
center. This galaxy has X-ray detection by the CXO with a 
HR of -0.55, and radio detection with /y^ =79 yuJy at 1 .4Ghz. 
A Spitzer MIPS flux density of 1580yuJy at 24yum gives it 
an LjR=7 x 10^ ^L© putting it in the category of Luminous 
Infra Red Galaxies (LIRGs). All the evidence points to the 
presence of an AGN. 

907047 This galaxy is an ongoing merger and morphologi- 
cal classification was not possible. At least two nuclei are 
clearly seen in the summed image, and the colormap shows 
the galaxy to be rather blue. 

907305 These are actually two merging galaxies with a sep- 
aration of about 0.35 arcsec. The fit failed. There are two 
distinct nuclei plus a halo that connects the two galaxies. 
The nuclei are bluer than the halo. 

907361 This is an early type bulge dominated galaxy. The 
galaxy is extremely red with a rest frame (U - B)o = 1.65 
mag. The galaxy disappears in the B band. 

907794 This galaxy is compact, isolated and featureless. The 
fit quality is very good and the appearance is that of an 
E/SO galaxy with red color as shown by the color map. The 
inferred rest frame (U - B)o = 1 .34 magnitude. This galaxy 
has a Spitzer MIPS /y^ = 1 04 Jy at 24yum yielding a L//? =2 x 
lO^^Lo making it a LIRG. The SFR for this galaxy is 35 
Mo/year. 

908243 A linear/edge on galaxy with a bright nucleus quite 
off'set from the center of the galaxy. The nucleus and the 
northern part of the galaxy are bluer than the rest of the 
galaxy. This is designated a tadpole. 

909015 This is a point source with the morphology undeter- 
mined. The object is red as seen from the color map. The 
integrated B-z color is ~5 magnitudes. Also, with a Mb of 
only -20.25, it is unlikely to be a quasar. Since this is a non- 
detection in Spitzer MIPS, it has an upper limit of 80 yuJy at 
24yum. Perhaps this is a star misidentified by VVDS. 

909093 This galaxy is shaped like a tadpole. There is an inten- 
sity maximum located towards the head of the tadpole (due 
south-west), ie not at the center of the galaxy. The fit is un- 
reliable due to the peculiar morphology. Again, the head of 
the tadpole is bluer than the rest of the galaxy. 

909429 This is a possible ongoing merger. At least two nuclei 
are distinctly visible in the i-band image. The nucleus at 
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the northern part of the system is blue, whereas the other 
nucleus is red. The fit failed. 

910413 This galaxy is shaped like an elongated tadpole. Two 
distinct intensity maxima are clearly seen. Both the maxima 
are bluer than the rest of the galaxy. The fit is unreliable. 
This might be an ongoing merger M2. This galaxy has a 
Spitzer MIPS /vo=248yuJy at 24yum, Lm = 1.6 x lO^^L©, 
making it a LIRG with a SFR of 27M0/year. 

911747 This is a possible ongoing merger system M2. Two 
intensity maxima are distinctly seen in the z band image, 
embedded in an elongated structure(visible in the summed 
image), which might be the tidal tail. Interestingly, one of 
the maxima is blue while the other is red. The fit failed. 

911780 This looks like an edge on disk dominated galaxy. The 
B band image looks like a tadpole and is markedly off'set 
with respect to the z band image. This is obvious from the 
color map, which shows that the blue region is off'set from 
the center of the galaxy. This might be a result of dust at- 
tenuation in the B band or alternatively off' center star for- 
mation. 

911843 A double nucleus is suspected in this system. The 
galaxy is shaped like a kidney. The north-west part of the 
galaxy is brighter and bluer than the rest of the galaxy. This 
brighter part is in fact composed of two nuclei which are 
barely resolved and are extremely blue with the rest frame 
(U - B)o = 0.49 magnitude. In some pixels the B-z color is 
negative, but this could partly be a consequence of the ACS 
z band red leak which transfers z band ffux away from the 
center of bright nuclei and makes the center look bluer than 
it really is (see Appendix B for details). This is a possible 
merger system. The fit failed. 

912744 A face on disk galaxy with spiral arms clearly dis- 
cernible in the summed image. There is a very bright and 
blue HII region on the south-west side. The patchy light 
distribution due to the presence of HII regions as well as 
the spiral arms make it difficult to obtain a meaningful fft. 
The colormap clearly shows the HII region to be bluer than 
the rest of the galaxy. 

913482 This is a compact, isolated and featureless E/SO galaxy 
with an AGN at the center. The core is point like and bluer 
than the outer regions of the galaxy. The fft is unreliable 
up due to the presence of the central point source. A vi- 
sual examination of the VVDS spectra of this galaxy con- 
ffrms that this is a type 2 AGN. This galaxy is a X-ray 
source with HR=0.52, as well as a Spitzer MIPS source 
with = 412.0 yuJy at 24yum and Ljr = 2.7 x lO^^L© mak- 
ing it a LIRG. Since this is an AGN, our prescription for 
calculating the SFR only gives an upper limit of 46Mo/year. 

914038 The morphology of this galaxy is knotty and it is prob- 
ably a merger system M2. Double nuclei is suspected at the 
south-east corner of the galaxy evident from the i band im- 
age. The color map shows a blue region on the south-west 
side of the galaxy. This galaxy is a Spitzer MIPS source 
with = ISO.OyuJy at 24yum and Ljr = 1.2 x lO^^L© mak- 
ing it a LIRG with a SFR of 20Mo/year. 

914895 This is most probably an ongoing merger. There is 
possibly an underlying disk galaxy with another galaxy 
merging with it (due north). The nucleus of the secondary 



member is very blue as seen in the colormap as well as the 
B band image. 

915400 This is a ffocculent spiral galaxy seen edge on as con- 
ffrmed by the fft. The color map shows blue colors on the 
outer edges of the disk, which might be indicative of inside- 
out disk formation. 

916137 This is a linear and ffocculent galaxy. The light dis- 
tribution is patchy and the center is not obvious. The fft to 
this irregular galaxy is unreliable. The north-east side of the 
galaxy is bluer compared to the other parts of the galaxy as 
seen in the color map. The summed image shows a more 
robust structure of the galaxy. 

916446 This system is a possible ongoing merger. There is a 
strong peak at the south-east side and a very faint peak at 
north-west side as seen in the i band image. There is some 
saturation in the pixels in the B band image, hence we have 
to disregard the color map for this galaxy. The fft failed. 

916866 This is an intermediate type galaxy (SO). The nucleus 
is slightly off'set from the center of the galaxy. The col- 
ormap shows that the south-west side is bluer than the rest 
of the galaxy. This galaxy is the host galaxy for the type la 
SN 2002ga (Strolger et al. [20Q4b . 

918147 This system is an ongoing merger Ml. The two nu- 
clei are distinctly visible in both B band and z band. The fft 
failed. The colormap shows that the south-east side nucleus 
is blue, and the north-west side nucleus is redder. This sys- 
tem has a Spitzer MIPS detection with /v^ = 206.0yuJy at 
24yum giving it an Lir = 3.9 x 10^ ^L©, making it a LIRG 
with a SFR of 67Mo/year. 

919573 This is an isolated galaxy with a suspected AGN due to 
the presence of a bright point like nucleus. The fft suggests 
this to be an SO galaxy, but is unreliable due to suspected 
AGN activity. The colormap shows that the south side of 
the galaxy is redder. This might be indicative of off'-center 
star formation in the north side of this galaxy. The inten- 
sity maximum in the z band is off'set from the center of the 
galaxy. This galaxy is an CXO X-ray source with HR=- 
0.40, and has a Spitzer MIPS detection with fy^ = 88.8 yuJy 
at 24yum, giving it an Ljr = 5.9 x lO^^L©. Due to the sus- 
pected AGN activity, our prescription for calculating the 
SFR gives an upper limit of lOMo/year. 

919595 This seems to be a disk galaxy with evidence of patchy 
star formation as seen in the color map as well as in the B 
band image. At least four prominent star formation HII re- 
gions can be seen in the B band image. The fft is unreliable 
due to the patchy nature of this galaxy. The z band image 
shows a single nucleus, hence the data does not suggest a 
recent merger. 

920435 Galfit suggests this to be an SO. The color map shows 
double nuclei (both are blue). The double nucleus is also 
seen in the B band image, but in the z band, the second 
nucleus disappears. This might be due to the presence of a 
dust patch which causes attenuation in the B band but not 
in the z band. Otherwise the galaxy is isolated with no signs 
of interaction/merger. Designated an SO. 

921406 The fft is excellent and this is designated to be a disk 
galaxy. It is an edge on disk galaxy with a low surface 
brightness extension whose reason is not clear. The whole 
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galaxy is blue in color. The disky nature of the galaxy is 
most clear in the summed image. 
922675 The galaxy is elongated. There is only one nucleus 
with a light distribution that is somewhat elongated to one 
side of the nucleus. The "tail" is bluer than the nucleus of 
the galaxy. The fit is excellent and the galaxy is designated 
to be a SO. 

922733 There is probable evidence for two nuclei specially in 
the summed image. The fit failed. This object is shaped like 
an elongated chain with multiple (two) peaks in intensity. 
This is possibly a merging system. The whole chain is blue 
in color. 

922761 The z band image of this galaxy is irregular with no 
obvious center. The fit is unreliable. The color map is blue 
throughout but is patchy with no obvious central intensity 
maximum. This is a blue irregular galaxy. 

923085 An irregular galaxy with multiple intensity maxima 
seen specially in the summed image. The fit is unreliable. 
There is a red patch at the core of the galaxy(dust?) sur- 
rounded by a relatively blue region. This galaxy has a 
Spitzer MIPS detection with /y^ = 84.4 yuJy at 24yum giv- 
ing it an LiR = 1.6 x lO^^L© making it a LIRG with a SFR 
of 26Mo/year. 

923926 Very compact, isolated and featureless galaxy. This 
galaxy has CXO X-ray detection with a HR of 0.12 which 
might be indicative of AGN activity, even though there are 
no discernible signs of point like activity from visual ex- 
amination of the galaxy image. It is designated SO but the 
fit is unreliable due to the possible presence of the AGN 
activity. The galaxy as a whole is red with a rest frame 
(U-B)o = 1.22 mag. 

924881 This galaxy has extensions at low surface brightness 
levels. The nucleus is off'set from the center of the galaxy. 
The fit is reasonable and this is designated to be an E type 
galaxy. The colormap shows the center of the galaxy to be 
bluer compared to the rest of the galaxy. 

926109 A simple visual examination of the B band image con- 
firms that this is an ongoing merger. There are at least four 
intensity maxima (star forming regions) with no tendency 
to form any regular pattern. These star forming regions also 
show up in the color map. The fit failed. 

926217 This is an isolated lenticular galaxy. The fit is excel- 
lent. The core of the galaxy is redder than the rest of the 
galaxy. There is some indication of a faint outer disk in the 
summed image, but that is not confirmed. This is designated 
as a lenticular. 

6. Number density evolution 

As has been widely reported, LCGs constitute a very rapidly 
evolving galaxy population in the intermediate redshift range 
(Phillips et al. fTWTl Werk et al. [2QQ4b . In order to verify this, 
we derived a rough estimate of the comoving number density of 
LCGs in the redshift range 0.5 < z < 1.2. After correcting for 
spectroscopic incompleteness using published values of Target 
Sampling Rate (TSR) by Ilbert et al. (2004), we get a value 
of 1.86 X 10"^ Mpc"^. Assuming a Poisson distribution yields 
a Icr error of -15% for our sample of 39 galaxies. We must 



be careful in comparing this number density with other pub- 
lished work because of the diff'erent sample selection criterion 
employed by diff'erent groups for choosing LCGs, especially in 
the absolute magnitude cut off'. With this caveat in mind, this 
number is comparable to the value reported by Phillips et al. 
([T997 ) as 1.2 xlO"^ Mpc"^ in the redshift range 0.7 < z < 1.0 
with an absolute magnitude cutoff' of Mb < -20. 

In order to look for evolution, we must compare our number 
density with the known number density of such sources in the 
local universe. Werk et al. ( 2004) have published the number 
density of local luminous compact blue galaxies (LCBGs) to be 
3.2 X 10""^ Mpc"^ for Hq = 70 km sec"^ Mpc"^ and an absolute 
magnitude cutoff' of Mb < - 19. 12. If we restrict their sample to 
only objects with Mb < -20 to agree with our sample selection 
criterion, we get a number density of 1 .4x10""^ Mpc"^, with the 
caveat that the Werk et al. sample suff'ers from small number 
statistics (only 5 objects) at such an absolute magnitude cutoff'. 

When we similarly restrict our sample to rest frame B - 
V < 0.6 to yield only blue LCGs in the redshift range 0.5 < 
z < 1.2, we are left with 33 of our 39 objects. This yields a 
comoving number density of 1.57 x 10"^ Mpc"^. A comparison 
with the number density of LCBGs in the local universe yields 
a factor ~ 1 1 fall in the comoving number density for blue LCGs 
from intermediate redshifts to the current epoch. The Poisson 
error on this number will be around ~ 50% which is mainly 
dominated by the small sample size at redshift zero. With this 
caveat of large errorbars, we still get a constraint on the number 
density evolution of blue LCGs which is consistent with earlier 
published results. 

7. Discussion 

As reported by earlier workers in this field (Guzman et al. ll9971 
Hammer et al. l2001 l). by their very nature LCGs are so compact 
that they yield little morphological information at intermediate 
redshifts, even using HST/WFPC2 images. Some noteworthy 
exceptions are the HST followup observations of CNELGs by 
Koo et al. (119941) , who found evidence for exponential light 
profiles for these objects. Also, Guzman et al. (1998) used 
HST/WFPC2 images of five CNELGs to reveal the presence 
of blue high surface brightness knots surrounded by a diff'use 
"exponential like" component. The knots were identified as the 
location of the current star formation and the diff'use exponen- 
tial like component was interpreted to be an older underlying 
population. 

In this work, we have performed quantitative morphologi- 
cal analysis of a sample of 39 LCGs using full 2 dimensional 
surface brightness profile fitting of the galaxy images. The ben- 
efit that we have, which was not available to our predeces- 
sors, is the technological advancement off'ered in the shape of 
ACS on board the HST. One of the advantages off'ered by the 
ACS imaging is the improved pixel scale of 0.05 arcsec/pixel 
which is a factor of ~2 better than the pixel scale of WFPC2, 
while at the same time having a larger footprint on the sky, 
which makes it ideal for large surveys. This better sampling be- 
comes specially important for compact sources such as LCGs 
whose typical ri/2 < 0.5arcsec, in which case the whole ex- 
tended emission (eg. disk) may be sampled by just a few pix- 
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els. This advantage of ACS has been fully leveraged by the 
HST/ACS GOODS survey which has the unique combination 
of depth, finer sampling of 0.03 arcsec/pixel and large area cov- 
erage which yields a large enough sample of LCGs to do mean- 
ingful statistical analysis. There have been other surveys which 
have some (but not all) of these three qualities, e.g. the Hubble 
Deep Field (WilHams et al. 1 1996b was deeper than GOODS but 
with an area coverage of only -5% of GOODS. GEMS (Rix et 
al. 12004 ) has larger area coverage than GOODS and the same 
pixel scale but is rather shallow (depth of single orbit). The 
HST/ACS Ultra Deep Field (Beckwith et al. l2QQ6l is much 
deeper than GOODS, but with area coverage which is much 
smaller (although see Noeske et al. (120061) , who derive a rea- 
sonable sample size of LCBGs by using photometric redshifts 
and a less stringent absolute magnitude cutoff' Mb < -18.5). 
Using the ratio of UDF area to GOODS area and our sample se- 
lection criteria for selecting LCGs, we would have got a sample 
of ~3 LCGs in the HST/ACS UDF. 

In this study, we find that LCGs constitute a significant 
-26% of the Mb < -20 galaxy population in the redshift range 
0.5 < z < 1.2. We derive SFR values ranging from a few to 
~ 65 Mo/year as expected for this class of objects. We estab- 
lish that our sample LCGs are intermediate mass objects with 
stellar mass ranging from 9.44 < logio(M/Mo) < 10.96, with a 
median mass oflogi^iMIMo) = 10.32. We also obtain a factor 
-11 fall in the comoving number density of blue LCGs from 
this redshift range to the current epoch, but this number is sub- 
ject to large uncertainities given the small sample size at zero 
redshift available from the literature. This implies that LCGs 
constitute one of the most rapidly evolving galaxy populations 
at intermediate redshifts. A plausible reason for the rapid de- 
cline in the number densities seems to be the extinguishing of 
the star formation and the associated fading, as pointed out by 
earlier workers (see Koo et al. ll995 , Guzman et al. ll998b . The 
eventual quenching of the starburst might result from super- 
novae driven galactic winds which might remove any residual 
gas from the galaxy (Guzman et al. 1 1996b . However, the nature 
of the descendants of these intermediate redshift LCGs will de- 
pend in a major way on the timing of this quenching and the 
possibility of subsequent bursts of star formation. It has been 
proposed (Koo et al. 1995, Guzman et al. lT996t that CNELGs 
should fade within a few Gyr. by as much as 4-7 magnitudes 
in the rest frame Mb to resemble today's spheroidals. However, 
one must be careful in comparing the CNELGs with our LCGs. 
A simple look at the stellar masses derived by us for our sam- 
ple of LCGs (Table [5]), make it clear why such a comparision 
is incompatible. CNELGs are mostly low mass systems, with 
M ~ 10^ Mo (Guzman et al.'1996'), whereas our LCGs are quite 
massive systems, with some of them having masses compara- 
ble to that of the Milky Way already in place by intermediate 
redshifts. Since they already have masses well in excess of that 
expected for dwarf galaxies, we can safely rule out the possibil- 
ity of these LCGs evolving into dwarf spheroidals in the local 
universe. Given the stellar masses that these objects have at in- 
termediate redshift, and the high star formation rates that we 
derive for these objects, they are clearly progenitors of inter- 
mediate mass objects in the local universe, out of which -70% 
are spirals in the current epoch (Hammer et al. 12005b . 



We find a small but significant ~22% of the LCGs to be disk 
dominated systems. The median disk scale length for the disk 
dominated LCGs is found to be -2.3 kpc. Interestingly, the col- 
ormaps for some of these disk dominated LCGs (eg. 905983, 
912744, 915400, 919595) show blue colors on the outer re- 
gions of the disk which might be indicative of inside-out disk 
formation as suggested by Hammer et al. (|2005b . However, 
equally interesting is the fact that none of these galaxies have 
a detection in Spitzer. This still allows for a SFR of several 
Mo/year (see Table [5]). However, the final fate of these LCGs 
will depend on when the star formation is quenched. Given the 
stellar masses that these objects have already assembled by in- 
termediate redshifts, these disk dominated LCGs are evidently 
the progenitors of intermediate mass disk galaxies in the lo- 
cal universe, as has been previously suggested by Phillips et al. 
([T997 I. 

Another relevant input which will have a strong inffuence 
on the final form of the descendants of these intermediate red- 
shift LCGs, is the evolution with redshift of the merger rate 
of galaxies, regarding which there has been wide debate in the 
community. Parameterizing the fraction of galaxies existing in 
merger systems at any given redshift as f(z) = /(0)(1 -h z)^, the 
reported values of a vary from < a < 4 (Le Fevre et al. 120001 
Patton et al. .2002. Consefice et al. [20031 Lin et al. (20041 Bundy 
et al. [20041 Bell et al. [20061 Lotz et al. [2006b . Even though 
the value of the slope a is still debatable, the merger rate es- 
timates are in reasonable agreement at least around redshift 
0.5-0.6 (see Bell et al. 2006). A large value of a would imply 
that the merger rate of galaxies was higher at higher redshifts. 
We have performed a calculation, along the lines of that per- 
formed by Le Fevre et al. ([2QQQb . using the HST/ACS GOODS 
dataset and found a rapid evolution in pair fraction of galaxies 
with redshift out to z~1.2 (Rawat et al. [in preparation| ). Even 
though this debate is far from settled, our work seems to in- 
dicate that a large fraction of L* galaxies in the local universe 
must have undergone a major merger, possibly triggering a new 
burst of star formation within the last 8 Gyr. In the present work 
we find that a large fraction (-36%) of the LCGs are in merg- 
ing systems, even though only about l/3rd of these are obvious 
mergers (Ml). These mergers have been identified essentially 
visually by the presence of multiple maxima of comparable in- 
tensity in the rest frame B band images. This was further aided 
by the use of color maps to distinguish multiple nuclei (gen- 
erally redder) from HII regions (generally bluer). Our defini- 
tion of merger is not directly comparable to those reported by 
other workers such as Lotz et al. (12006b who use a quantita- 
tive definition to pick up merging systems. Despite this caveat, 
our obvious merger fraction (Ml) (~ 12%) is quite compara- 
ble to the typical merger fraction reported for the general field 
galaxies by Lotz et al. (12006b , who find the fraction of galaxies 
with Mb < -20.5 that are in merging systems to be -7% at 
intermediate redshifts. If we also include our possible merger 
systems, this is evidence in support of an unusually high merger 
rate for LCGs, and points towards a dynamic evolutionary his- 
tory of LCGs in the intermediate redshift range. It is plausible 
that LCGs might represent a transient, starbursting (possibly 
merging) phase in the evolution of a galaxy. In that case, the 
observed decline in the number density of LCGs from z~ 1 to 
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the current epoch can be naturally explained in terms of the 
decline in the merger rate of galaxies over the same redshift 



range (Rawat et al. in preparation). This is in agreement with 
the scenario suggested by Hammer et al. (^005 ) in which the 
LCGs are just the merger phase in the hierarchical growth of a 
galaxy. 

The eventual fate of the LCGs undergoing merger is still 
an open issue. Hammer et al. (l2QQ5t had suggested that there 
must be a mechanism at work that rebuilds the disk in a galaxy 
after a major merger event, perhaps fueled by inf ailing gas left 
over after the merger. They had argued that if the fate of every 
LCG is to end up as an early type galaxy, it will overproduce 
the number of early types in the local universe (see Table 2 
in Hammer et al. l2QQ5l) . In this scenario, we can envisage that 
some of the LCGs that are seen in merger phase will go on to 
form a disk as they evolve to z=0. Recent work on numerical 
simulations of galaxy mergers (Robertson et al. l2QQ5b indicates 
that gas rich mergers at high redshift can lead to the formation 
of rotationally supported disks in merger remnants if feedback 
mechanisms limit the conversion of gas into stars. This gas rich 
merger driven disk galaxy formation at high redshifts degener- 
ates into merger driven elliptical galaxy formation at lower red- 
shifts, as star formation depletes the available gas and makes 
gas rich mergers less likely. 

Another support for this scenario comes from the new strin- 
gent constraints being put on the evolution of early type galax- 
ies (ETGs) since z ~ 1 by Cimatti et al. ('2006). They have 
proposed that at every redshift, there is a critical luminosity (or 
mass) threshold above which virtually all ETGs seem to be al- 
ready in place. This is also known as the downsizing scenario. 
Now, given the fact that our sample LCGs are intermediate 
mass objects with a mass range 9.44 < Logio(M/MQ) < 10.96 
(see Table [5]), this does not seem to favor their evolution into 
early types at the current epoch, simply because most of the 
ETGs in this mass range were already in place latest by inter- 
mediate redshifts (see Fig. 2 in Cimatti et al. 2006) . 

Interestingly, Lotz et al. (120061 ) have claimed a large fac- 
tor -2.5 increase in the number density of early types E/SO/Sa 
from redshift z~l.l to z~0.3 for a volume limited sample of 
Mb < -20.5. However, a simple look at their Table 2 shows 
that most of their evolution comes from the highest redshift bin, 
with the lower redshift data consistent with little or no evolu- 
tion out to a redshift of ~ 1 . Their result is in fact similar to the 
result published by Smith et al. (I2005L who find the fraction of 
early type galaxies to remain constant with Z^+^o = 0-4 ±0.1 
till z ~ 1 in the low density regions. This is also consistent 
with Cimatti et al. (12006b result that most of the massive (lu- 
minous) ETGs are already in place by intermediate redshifts. 
Now, given the fact that LCGs constitute only ~ 1/4 of the field 
galaxy population with Mb < -20, it is entirely possible that a 
large fraction of LCGs might still evolve into ETGs by the cur- 
rent epoch, without upsetting the constant Z^+^o fraction con- 
straint discussed above within their quoted errorbars. So, even 
though we cannot rule out an evolutionary scenario of inter- 
mediate redshift LCGs evolving into early types in the current 
epoch (given the uncertainties of the early type fraction evolu- 
tion with redshift), the current results on ETG fraction evolu- 
tion do not seem to encourage a large scale evolution into ETGs 



in this luminosity range, over the redshift range under consid- 
eration. This result will become stronger, as the errorbars on 
the Z^+^o evolution with redshift improve with larger surveys 
becoming available. 

If we combine all the evidence listed above, we envisage 
that it is plausible that some of the LCGs that are undergo- 
ing mergers at intermediate redshifts might go on to form a 
disk from gas left over from the merger event, by the time they 
evolve to the current epoch. 

Puech et al. (12006b have recently reported on the morpho- 
logical mix of a sample of 17 intermediate redshift LCGs using 
the technique of mapping the velocity field of a galaxy with the 
Integral Field Unit (IFU) spectrograph FLAMES/GIRAFFE on 
the VLT. They report a figure of -18% galaxies being consis- 
tent with rotating disks compared to the -22% in this work. 
Furthermore, they report that -82% of the LCGs in their sam- 
ple have perturbed or complex kinematics, with roughly half 
of them being probable mergers, a number which again agrees 
with our result designating -36% of the LCGs as merging sys- 
tems. The techniques applied by Puech et al. are entirely diff'er- 
ent from ours: 

a) They probe ionised gas, whereas we probe stellar mass, 

b) They use ground based spectroscopy while we use space 
based imaging. 

The fact that we get results that are comparable gives us in- 
creased confidence in our results. 

Noeske et al. ( 2006 ) have recently reported on the morphol- 
ogy of LCGs from the Hubble Ultra Deep Field in the redshift 
range 0.2 to 1.3. They find typical scale lengths of the extended 
components to be < 2 kpc, and conclude that these objects will 
evolve into small disk galaxies and low mass spheroidal or ir- 
regular galaxies in the local universe. However, their sample se- 
lection criteria are not fully compatible with those of Hammer 
et al. (12001 ), which we have used in this work, and which might 
favor progenitors of larger local galaxies. We find that the me- 
dian disk scale length for the disk dominated LCGs in our work 
is -2.3 kpc (see Table [5]). These LCGs are intermediate mass 
objects and given the limits on their SFR of a few Mo/year, are 
likely to be the progenitors of intermediate mass disk galaxies 
in the local universe. 

8. Conclusions 

An eff'ort was made to obtain quantitative morphological pa- 
rameters for a representative sample of 39 intermediate redshift 
(0.5 < z< 1.2) luminous compact galaxies (LCGs) and to clas- 
sify them based on their rest frame B band morphology and B-z 
color maps using multiband imaging data from the HST/ACS 
GOODS survey. The main result that we obtain for their mor- 
phological mix is as follows: Mergers: -36%, Disk dominated: 
-22%, SO: -20%, Early types: -7%, Irr/tadpole: -15%. 

We find that LCGs constitute -26% of the Mb < -20 
galaxy population in the redshift range 0.5 < z < 1.2. We 
also obtain a factor - 1 1 fall in the comoving number density of 
blue LCGs from the redshift range 0.5 < z < 1.2 to the current 
epoch, using published values for the comoving number den- 
sity of blue LCGs in the local universe (Werk et al. .2004 ). We 
establish that our sample LCGs are intermediate mass objects 
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with stellar mass ranging from 9.44 < logiQCM/M©) < 10.96, 
with a median mass of log^QCM/M©) = 10.32. We also derive 
SFR values ranging from a few to ~ 65 Mo/year as expected 
for this class of objects. 

We find that the median disk scale length for the disk dom- 
inated LCGs in our work is ~2.3 kpc. Some of the disk domi- 
nated LCGs show blue colors on the outer regions of the disks 
which might be indicative of inside-out disk formation as sug- 
gested by Hammer et al. ( 2005i) . These disk dominated LCGs 
are intermediate mass objects, and given the star formation 
rates that we derive for these objects, they are clearly progeni- 
tors of intermediate mass disk galaxies in the local universe. 

We also establish that we cannot rule out an evolutionary 
scenario comprising of intermediate redshift LCGs evolving 
into early types in the current epoch, given the typical erro- 
bars on the early type fraction evolution available from litera- 
ture. Despite this, we note that the current results on ETG frac- 
tion evolution, do not seem to encourage a large scale evolution 
into ETGs in this luminosity range, over the redshift range un- 
der consideration. We envisage that some of the LCGs that are 
classified as merging systems, might go on to rebuild their disks 
and evolve into disk galaxies in the local universe. 
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Appendix A: Determination of absolute magnitudes 

We have used the technique of fitting template galaxy spectrum 
to a set of photometric datapoints for each galaxy to obtain the 
spectral energy distribution (SED) of the object, which can then 
be used for deriving the K-correction. Obviously, as in the case 
of any fitting procedure, the greater the number of photometric 
datapoints, the better we will be able to constrain the galaxy 
SED. Hence, in addition to using the four band HST/ACS pho- 
tometric catalog for the GOODS/CDFS sources, we obtain two 
more photometric datapoints using the J & Kg imaging data 
from the ESO GOODS/EIS release. SExtractor based photom- 
etry for the J & Ks fields was performed by us using the photo- 
metric zeropoints provided as part of the GOODS/EIS release: 
vl.O. The photometric catalogs for the J and Kg band in this 
field have -8000 sources each. The diff'erential number counts 
for sources in the Kg band is given in Fig.[9l 

In the absence of detailed completeness curves for the Kg 
band source catalog, the diff'erential number counts can be used 
as a useful indicator for the completeness of the Kg band cata- 
log. As we can see from Fig. [9] the number counts are steadily 
rising till Ks -24 magnitude, and it starts to taper off' for fainter 
magnitudes. Hence we can assume that the Ks band catalog 
suff'ers from incompleteness for sources fainter than Ks ~24 
magnitude. Since in our analysis work, we are only interested 
in galaxies having spectroscopic redshifts from the VVDS sur- 
vey, all of which have Iab < 24, we will not miss any object in 
the Ks band, unless it has an extremely unusually blue color of 
(I-Ks)<0. 



log N(m) Vs m plot for GOODS/EIS Ks band 




Having compiled the list of all the objects in the CDFS hav- 
ing photometry in 6 bands plus spectroscopic redshifts from 
the VVDS as mentioned above, we derived the absolute mag- 
nitudes for those objects as follows. 

The transformation equation relating observed frame pho- 
tometry with rest frame absolute magnitudes is... 



mB = Mb + DM(z) + kB(z) + Ab 



(8) 



18 20 22 

Ks magnitude 



Fig. 9. Ks band number counts for the VLT GOODS/EIS 



where, 

mB : apparent magnitude in B band 

Mb : absolute magnitude in B band 

DM(z) : distance modulus 

kB(z) the K-correction in B band 

Ab : Galactic extinction correction in B band 

For the purpose of fitting template galaxy spectra, we used 
a code developed by Hammer et al. (120011) . To derive ab- 
solute magnitude at a given rest-frame wavelength Arest^ we 
consider two observed frame wavelengths Aobsi & ^obsi^ such 
that ^^obsiliX + z) & ^obsiliX + z) encloses the rest-frame Xrest^ 
where z is the redshift of the object. To this set of galaxy ob- 
served frame magnitudes, we fit a galaxy spectral template de- 
rived from stellar population synthesis models (Bruzual and 
Chariot 1199312003ft . For our fitting purpose, we used an ex- 
ponentially decreasing SFR model with e-folding timescale of 
T = 0.5 Gyr. We then vary the age of this model to optimize 
the fit to observed frame magnitudes. The template has been 
optimized to minimize the residuals between the actual galaxy 
magnitudes and the galaxy magnitudes reconstructed from the 
galaxy SED fit. We then use an interpolation to calculate the 
ffux at the given rest- frame Xrest- This code then uses a specified 
cosmology to calculate the distance modulus and output the ab- 
solute magnitudes directly. We used a fiat ACDM cosmology 
with n^=0.3 and Qa=0.7 and = TOkmsec"^ Mpc"^. 

We used the Galactic extinction maps published by 
Schlegel et al. (119981) to estimate the Galactic extinction in 
each of the four HST/ACS photometric bands for the CDFS. 
Since the CDFS is a high Galactic latitude field, we found 
that the extinction values were quite small; being of the order 
of photometric accuracy of our catalog i.e. -1-2%. Hence 
we decided not to use these extinction values for our final 
estimation of the absolute magnitudes. 

Appendix B: Red Halo in the HST/ACS F850LP PSF 

It has been well documented and reported that at redder wave- 
lengths, all HST/ACS CCD detectors suff'er from an extended 
halo due to scattered light (Sirianni et al. 12005 , Gilliland and 
Riess 2002 ). This eff'ect is particularly severe in the F850LP fil- 
ter and depends upon the color of the PSF star being observed. 
This happens due to the fact that the CCD becomes more trans- 
parent to photons of longer wavelength. This means that for 
a redder star, the scattering of light into the halo is more se- 
vere compared to a bluer star. This variation of the PSF with 
the color of the PSF star can be quite large and can potentially 
aff'ect the structural parameters that we derive for our sample 
galaxies. 

In order to check for this possibility, we plotted the F850LP 
filter half light radius R1/2 of the PSF stars used by us against 
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the B-z color of the PSF star. This is shown in Fig.[TOl Different 
symbols have been used for faint (m^ > 18.5) and bright (m^ < 
18.5) stars using an arbitrary definition of what constitutes a 
bright star. 



3-z color Vs r^/2(zband) of 64 PSF stars m CDFS 



m >18.5 



^ o 



2 4 
-z color of PSF star 



Fig. 10. Variation of half light radius R1/2 of the PSF stars (in 
the F850LP filter) used by us against the B-z color of the PSF 
stars. Diff'erent symbols have been used for faint (m^ > 18.5) 
and bright (m^ < 18.5) stars using an arbitrary definition of 
what constitutes a bright star. Notice that at a given color, the 
brighter stars tend to have larger R1/2 (even though there is a 
considerable amount of scatter). 



Notice in Fig. [TOl that at a given color, the brighter stars 
tend to have larger R1/2. Also notice the small gradient towards 
increasing R1/2 for redder stars. However, it is seen that even 
though there is a general trend for redder stars to have slightly 
larger 7? 1/2, this eff'ect is < 5%. The only catastrophically large 
7? 1/2 is seen in relatively bright stars. From Fig.[TOlwe see that 
only four stars in our whole sample show this anomolous be- 
haviour. We also find that all of the "peaks" that we had seen 
in Fig. [3] belong to these same red and bright PSF stars. Once 
these four rogue stars were identified, there were removed from 
our list of PSF stars and were not used in any of our analysis. 
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Table 2. Galfit results in F850LP band for those LCGs which are well fit by bulge+disk two-dimensional structure. 



ID 


Bulge Mag^ 


Bulge re 
(pixels)^ 


Bulge n 


Disk mag^ 


Disk rd 
(pixels)^ 


^Red 




B/T 




904604 


23.66±0.02 


1.92±0.03 


1.21±0.08 


22, 


.87±0. 


,01 


10.56+0.15 


1 


.063 


0. 


,33+0. 


,01 


905632 


24.04±0.08 


2.56±0.31 


4.29±0.43 


23, 


.29±0. 


,03 


4.87+0.07 


1 


.093 


0. 


,33+0. 


,03 


905983 


24.55±0.23 


19.52±12.92 


12.01+3.70 


21, 


.92±0. 


,00 


8.00+0.06 


1 


.129 


0. 


,08+0. 


,02 


907361 


21.63±0.01 


5.86±0.11 


2.71±0.04 


23, 


.42+0. 


,07 


17.33+1.13 


I 


.259 


0. 


,84+0. 


,07 


907794 


22.76±0.11 


2.67±0.52 


4.64±0.70 


23, 


.53±0. 


,19 


8.05+0.46 


1 


.099 


0. 


,67+0. 


,09 


908243 


24.07± 0.04 


3.65±0.11 


0.82±0.05 


22, 


.21±0. 


,01 


10.46+0.07 


1 


.168 


0. 


,15+0. 


,01 


909093 


23.10± 0.02 


8.16±0.07 


0.53±0.02 


23, 


.76±0. 


,04 


11.33+0.41 


1 


.093 


0. 


,65+0. 


,02 


911780 


22.72±0.02 


7.44±0.05 


0.73±0.01 


21 


.49±0 


.00 


11.27+0.06 


1 


.185 


0. 


,24+0. 


,00 


913482 


22.04±0.00 


1.61±0.01 


0.98±0.01 


21 


.83±0, 


.00 


6.64+0.05 


1 


.111 


0. 


,45+0, 


,00 


915400 


23.55±0.42 


51.35±73.45 


17.69±8.14 


22, 


.25±0, 


.01 


10.27+0.08 


1 


.260 


0. 


,23+0, 


,01 


916137 


23.29±0.03 


14.44±0.12 


0.20±0.01 


23 


.73±0, 


.05 


6.83+0.17 


1 


.094 





.60+0, 


,02 


916866 


23.75±0.05 


11.88±0.30 


0.71±0.03 


23, 


.39±0, 


.03 


5.76+0.09 


1 


.137 





.42+0. 


,03 


919573 


22.13±0.05 


8.68±0.90 


7.43±0.35 


22, 


.25±0, 


.03 


10.89+0.13 


1 


.131 





.53+0. 


,03 


919595 


























920435 


23.81±0.02 


3.36±0.03 


0.07±0.06 


23 


.77±0, 


.01 


4.38+0.14 


1 


.116 





.49+0. 


,01 


921406 


23.47±0.08 


7.03±0.55 


2.03±0.12 


22, 


.61±0 


.04 


9.90+0.12 


1 


.125 





.31+0. 


,03 


922675 


22.75±0.02 


17.14±0.27 


1.14±0.02 


22, 


.47±0, 


.01 


4.85+0.03 


1 


.089 





.44+0. 


,01 


922761 


24.33±0.05 


8.86±0.17 


0.23±0.03 


23, 


.05±0, 


.02 


7.64+0.10 


1 


.095 





.24+0. 


,01 


923085 


22.70±0.01 


8.25±0.08 


0.93±0.02 


24 


.75±0, 


.04 


3.51+0.14 


1 


.212 





.87+0. 


,01 


923926 


22.58±0.02 


3.31±0.19 


7.76±0.77 


23, 


.22±0, 


.05 


1.97+0.02 


1 


.113 





.64+0. 


,02 


924881 


22.59±0.01 


6.03±0.15 


2.99±0.07 


23, 


.93±0, 


.03 


2.07+0.05 


1 


.167 





.77+0. 


,01 


926217 


22.70±0.05 


11.09±1.03 


5.48±0.27 


22, 


.80±0, 


.03 


8.21+0.12 


1 


.097 





.52+0. 


,03 



^ Integrated Bulge/Disk magnitude given by galfit. 
^ I pixel - 0.03 arcsec 



Table 3. Galfit results in F775W band for those LCGs which are well fit by bulge+disk two-dimensional structure. The columns 
are the same as in Table [2l 



ID 


Bulge Mag 


Bulge 
(pixels) 


Bulge n 


Disk mag 


Disk rd 
(pixels) 


^Red 


B/T 


904604 


24.67+0.04 


1.83+0.07 


1.31+0.14 


23.91+0.01 


10.43+0.33 


1 


.088 


0.33+0.01 


905632 


24.24+0.06 


2.71+0.23 


3.29+0.27 


23.97+0.04 


4.36+0.09 


1 


.132 


0.44+0.03 


905983 


24.73+0.30 


19.08+17.34 


15.77+5.67 


22.19+0.00 


8.33+0.05 


1 


.200 


0.09+0.02 


907361 


22.26+0.03 


5.27+0.17 


2.50+0.05 


23.76+0.10 


14.45+1.43 


1 


.322 


0.80+0.04 


907794 


















908243 


24.45+0.03 


3.52+0.09 


1.01+0.05 


22.43+0.01 


10.85+0.07 


1 


.195 


0.13+0.00 


909093 


23.42+0.02 


11.87+0.23 


1.43+0.04 


23.89+0.02 


3.92+0.06 


1 


.195 


0.61+0.01 


911780 


22.81+0.01 


7.19+0.04 


0.76+0.01 


21.70+0.00 


11.16+0.05 


1 


.295 


0.26+0.00 


913482 


22.29+0.00 


1.74+0.01 


1.05+0.01 


22.22+0.00 


7.08+0.06 


1 


.221 


0.48+0.00 


915400 


25.00+0.34 


7.50+6.83 


11.23+4.06 


22.41+0.01 


11.03+0.07 


1 


.333 


0.08+0.01 


916137 


23.60+0.04 


13.92+0.10 


0.22+0.01 


24.27+0.07 


7.52+0.18 


1 


.124 


0.65+0.03 


916866 


















919573 


22.19+0.08 


41.48+10.85 


17.97+1.39 


22.68+0.02 


9.38+0.10 


1 


.172 


0.61+0.05 


919595 


24.51+0.01 


8.59+2.31 


0.05+0.04 


22.32+0.00 


7.66+0.04 


1 


.182 


0.12+0.00 


920435 


24.20+0.03 


4.99+1.02 


0.02+0.03 


24.17+0.02 


2.53+0.04 


1 


.143 


0.49+0.02 


921406 


24.53+0.12 


5.51+0.74 


2.63+0.29 


22.74+0.02 


9.68+0.09 


1 


.186 


0.16+0.02 


922675 


22.78+0.01 


17.64+0.22 


1.11+0.02 


22.66+0.01 


4.86+0.03 


1 


.142 


0.47+0.01 


922761 


24.81+0.03 


8.77+0.13 


0.13+0.03 


23.28+0.01 


8.38+0.09 


1 


.124 


0.20+0.01 


923085 


23.44+0.01 


10.05+0.08 


0.40+0.01 


25.66+0.04 


2.87+0.15 


1 


.248 


0.89+0.01 


923926 


23.26+0.04 


6.42+1.44 


10.47+1.89 


23.84+0.05 


1.79+0.03 


1 


.163 


0.63+0.03 


924881 


22.81+0.01 


5.94+0.17 


3.49+0.09 


24.18+0.03 


2.13+0.04 


1 


.277 


0.78+0.01 


926217 


23.00+0.09 


8.41 + 1.03 


3.64+0.23 


23.00+0.07 


10.65+0.25 


1 


.119 


0.50+0.05 
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Table 4. Galfit results in F606W band for those LCGs which are well fit by bulge+disk two-dimensional structure. The columns 
are the same as in Table [2] 



ID 


Bulge Ma.g 


Bulge Tg 
(pixels) 


Bulge n 


Disk iiia.g 


Disk 
(pixels) 


^Red 


B/T 


004604. 


Ifi ^9+0 OS 

z, vj . J) z, i vj . vj J 


I 99+0 12 


91+0 25 


7< S8+n 0^ 


1 7 77+0 81 


1 316 


^4+0 09 




74 Ql +0 07 


1.73+0.03 


71 +0 ns 


74 ^;q+0 01 


4 ss+n 14 


1.355 


4s+n ni 


















907361 
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Table 5. Catalog of derived parameters for LCGs 



Stellar Mass^ SFR(p) _Rest Frame B _ 
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SO 
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SO 
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^ Derived using a diet Salpeter IMF. In order to derive values using a classic Salpeter IMF, add 0.15 dex to the values listed here. 
^ X-ray hardness ratio, defined as (H-S)/(H+S) where H and S are the counts in the hard and soft bands respectively. 

Galaxy type — E: 0.8 < B/T < 1; SO: 0.4 < B/T < 0.8; S: 0.0 < BjT < 0.4; Irr: irregular; Tad: Tadpole. 
^ Q quahty factor — 1: secure; 2: possibly secure; 3: insecure; 4: fit failed. 
^ Exponential disk scale length(kpc) for disk dominated galaxies. 
^ Merging? — Ml: obvious merging; M2: possible merging. 

2 This object could not be fitted properly due to the presence of a strong point source at the center. 
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[. 8. and the B-z color map images. Explanation is 
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[. 8. Continued. 
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[. 8. Continued. 
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Fig. 8. Continued. FH50LP and the B-z color map images. For each target, the name and redshift are labeled at the top-left and 
top-right in the F^SOLP image. Our morphological classification is indicated at the bottom-left of each image. The size of each 
image is 30x30 kpc and the bar at the bottom-right of each image indicates 1 .0 arcsec. North is up, East is to the left. The intensity 
scaling is the same (squareroot) in all the grayscale images, while the dynamic range is diff'erent for diff'erent objects owing to 
the fact that the max or central intensity can be diff'erent for diff'erent objects. The color bar shows the B-z color range from to 
4 except for object ID 904604, 907361 and 909015 where it ranges from to 6. The blank in object ID 907305 is due to imaging 
at the edge of the GOODS CDFS region. 



